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Abstract. We implemented a quasi time-dependent 2D 
stochastic model of solar modulation describing the trans- 
port of cosmic rays (CR) in the heliosphere. Our code can 
modulate the Local Interstellar Spectrum (LIS) of a generic 
charged particle (light cosmic ions and electrons), calculating 
the spectrum at 1 AU. Several measurements of CR antiparti- 
cles have been performed. Here we focused our attention on 
the CR antiproton component and the antiproton over proton 
ratio. We show that our model, using the same heliospheric 



parameters for both particles, fit the observed ^ ratio. 



We 



show a good agreement with BESS-97 and PAMELA data 
and make a prediction for the AMS-02 experiment. 



1 Introduction 

Galactic cosmic rays (GCRs) are nuclei, with a small com- 
ponent of lep tons, mainly produced by supernova remnants 
I Blasil ( 120 10 )1. confined by the galactic magnetic field to 
form a isotropic flux inside the galaxy. Before reaching the 
Earth orbit they enter the heliosphere, the region where the 
interplanetary magnetic field is carried out by the solar wind 
(SW). In this enviroment they undergo diffusion, convection, 
magnetic drift and adiabatic energy loss, resulting in a reduc- 
tion of particles flux at low energy (<1-10 GeV) depending 
on solar activity and polarity. This effect is known as so- 
lar modulation. We have developed a 2 D (radius and heho - 
colatitude) model of GCR propagation IIBobik et al.l (l2003h l 
in the heliosphere, by using stochastic differential equations 
(SDEs). The model depends on measured values of the SW 
velocity on the ecliptic plane (Vq), tilt angle (a) of the neutral 
sheet and estimated values of the diffusion parameter (ko): 
details on parameters are discussed in section 2 and 3. This 
model includes drift transport due to magnetic field curva- 
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ture and gradients, as well the presence of a tilted neutral 
sheet describing properly periods of low and medium solar 
activity. Modulated fluxes depend on solar activity but also 



on par ticle charge and solar magnetic polarity [Boella et al 
(1200 ih l. 



2 Stochastic 2D Monte Carlo code 

The GCR transport in the Heliosphere is described by 
a Fokker-Planck equation, the so-called Parker equation 
[Parker (1965)]: 
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where U is the cosmic ray number density per unit interval of 
particle kinetic energy, t is the time, T is the kinetic energy 
(per nucleon), Vswi the SW velocity along the axis Xi, 
is the drift velocity related to the antisymmetric part of dif- 
fusion tensor [Jokipii and Levy (1977)], K^^ is the symmet- 
ric part of flie diffusi on tensor and g={T + 2To)/(T + Tq) 
iCl eeson et al. ( 1967 )]. where Tq is particle's rest ene rgy. 



This partial differential equation is equivalent IGardinen 



(119891) 1 to a set of ordinary SDEs that can be integrated with 
Monte Carlo (MC) techniques. The integration time step 
(At), is taken to be proportional to (r is the distance from 
the Sun) avoiding oversampli ng in the outer heliopshere an d 
therefore saving CPU time OAlanko-Huotari et al.l (120071) ] . 
We considered the 2D (radius and colatitude) approximation 
of Eq. [H and from this we calculate the equivalent set of 
SDEs: 
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where /x=cos0, with 9 colatitude, and Rg is a gaussian dis- 
tributed random number with unitary variance. Here the 
drift velocity is split in regular drift (radial drift vo,., 
tu dinal drift vna) and neutral sh eet drift (wdns) described 
by iHattingh and Burgen (119951) . Th e radial diffusion coeffi - 
cientis Krr=K\\cos'^il)+K x^rsin^ip IPotgieter. et al.l (1 19931) 1. 
where ifj is the angle between radial versor and direction 
of the solar magnetic field describe d below. The latitudi- 
nal co efficient is K0g=K±g [e.g., see iPotgieter and Le Roux 
d 19941) 1. We note as the perpendicular diffusion coef- 
ficient has two components, one in the radial direction 
(K I r) and one in the polar direction (K±g) as shown in 
Potgietei tOOd) . We define (if±)o. as the ratio between 
perpendicular and parallel diffusion coefficients, therefore 
K±^={K±)o,K\\. We fixed this value: (^4:^)0, =0.05 while 
{K±)ob = fl9){Kj^)o,. (where f{e) = 10 close to the p oles 
and f{6)^l in the equatorial regionl llPotgieted (l200ol) 1. to 
reproduce the correct magnitude and rigidity dependence 
of the l atitudinal cosm i c ray proton and e l ectron gradi- 
ents [cf. IPotgieter et all (Il997b : iBurger et al.l (l2000l) 1. The 
parallel diffusion coefficient is K^^=koPKp{P){B(^/3B) 
I Potgieter and Le Rouxl(ll994 l: here fcn - 0.05 - 0.3 x IQ-^ 



is a diffusion parameter depending on the 
solar actvitiy (see section |3]l, f3 is the particle velocity, P 
is the CR particle's rigidity, Kp=P, is the value of 
heliospheric magnetic field at the Earth orbit, and B is 
the magnitude of the Helio spheric Magnetic Field (HMF) 
I Hattingh and Burgeij d 1 995l) 1 : 



B: 
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where A is a coefficient that determines the field polarity and 
allows |B| to be equal to B(^, i.e., the value of IMF at the 
Earth orbit; 9' is the polar angle deter mining the position 
of the heliospheric current sheet (HCS) IJokip ii and ThomasI 
(Il98lh l: H is the Heaviside function, thus [1 - 2H{9 - 9')] 
for the change of sign between the two regions - above and 
below the HCS - of the heliosphere; finally F = tani/; ~ 
with tjj the spiral angle. The Parker field has 
been modified introducing a small latitudinal component 
Be=^{r/rQ)5{9) with 5{9) =^.7 y.lQ-^ / iiin9, thus allow- 
in g one to obtain V ■ B = and a field magnitude according 
to lJokipii and Kotal (Il989l) : 



(6) 



that increases the magnitude of the HMF in the polar regions 
without a modification of the field topology. This component 
produces a lower magnetic drift velocity in this region, with 
the effect of a lower CR penetratio n along polar field hnes in 



the inner part of the heliosphere llJokipii and Kota (Il989l) 1. 



We use a SW broad smoothed profile according to Ulysses 



data for periods of low solar activity I Mc Comas et al, 
(l2000l) 1. described by the relation Vsn, {d) = Vmax if 6^ < 30° 
or > 150° and V^^ (0) = • (1 + | cos6l| ) if 30° < 6* < 150° 
where Vb is approximately 400 km/s and Vmax is 760 km/s. 
Drift effects are included through analytical effective drift 
velocities: in the Parker spiral field we evaluated drift due to 
gradient, curvature and neutral sheet that modify the integra- 
tion path inside the heli osphere. We adopted t he app roach 
of Potgieter and Moraal llPotgieter and Moraall d 1985b ], be- 



cause it is able to reproduce the effects of drift in both quiet 
and activ e solar periods [a d iscussion on other models can be 



found in Bobiketal 



dlOlO")]. In this model the drift coeffi- 
cient is modified with a transition function that simulates the 
effect of a wavy neutral sheet. The sharpness of this func- 
tion is related to a angle, expanding or shrinking the region 
of influence of neutral sheet drift. As LI S, both for pro tons 
and antiprotons, we use the ones used in lCasaus (l2009l) and 
obtained from GalpropQ- 

3 Parameters and Data Sets 
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Fig. 1. fco values versus monthly SSN. The linear fit is also shown. 
Reported values cover the time period 1951-2004. 



Values of the tilt angle a are cornputed using two differ- 
ent models, described in lHoeksemal (119951) . fitting separately 



periods of increasi ng solar activity and period s of decreas- 
ing solar activity I Ferreira and Potgieterl (120041 )1. The three 
drift components do not depend on external parameters, ex- 
cept the solar polarity, so A>0 for positive periods and A<0 



1http://galprop.stanford.edu/webrun/| 



Davide Grandi: Antiproton modulation in the Heliosphere and AMS-02 antiproton over proton ratio prediction 



3 



for negative periods OBobik et al.l (120031) 1. We selected CR 
p and p data from several experiments in order to compare 
and tune model results. We modulated separately p and p 
LIS spectra and then we computed the ratio. In this paper 
we show experimental data taken during periods of low so- 
lar activity: the comparison with BESS-97 (A>0 July 1997, 



Orito et all [iQOO,)) and P AMELA (A<0 from 2007 to 



see 

2008, see Adriani et al. ( 2010l) ). Vq and values for these 
periods were obtained from NSSDC OMNIWeb systeniby 
27 daily averages , while tilt angle values from the Wilcox 
Solar Laboratory OHoeksemal (Il995h l . We estimated the val- 
ues of ko, needed to evaluate the CR modulation in differ- 
ent conditions, from the modulation parameter reported in 
Usoskin et al.l ( l2005b . We searched a relation between the es- 
timated ko values and the monthly Smoothed Sunspot Num- 
bers (SSN). We found that there is a nearly linear relation 
between ko and values (see Fig. [Hi, with a Gaussian 

distribution of the best fit with a RMS of 19%. This is a first 
crude estimation, we will perform a more complex anaysis, 
e.g. fitting separately different solar phases, in order to avoid 
systematics in the relation and to reduce the RMS. In this 
way we can use the estimated SSN values to obtain the dif- 
fusion coefficient feg. Following this approach we introduced 
in our code a gaussian random variation of fco with a RMS 
of 19%. Results of the simulation with and without the gaus- 
sian variation are consistent inside the indetermination of the 
code (around 5%). Our code simulates a diffusive propaga- 
tion of a CR entering the heliosphere from its oute r limit , 
that we located at 100 AU (note that in lOecker et al.l (l2005h 
the Termination Shock is located at 94 AU) , and reaching 
the Earth at 1 AU; the effects of heliosheath and termination 
shock are not taken into account in the present model. We 
evaluated the time tg^ needed by the SW to expand from the 
outer corona up to 100 AU, with a minimum speed of ^400 
km/s it takes nearly 14 months, while the time interval t^v of 
the stochastic evolution of a quasi particle inside the helio- 
sphere from 100 AU down to 1 AU is between 1 month (at 
200MeV) and few days (at 10 GeV). This scenario, where 
Tev < tsw and tsiu > > 1 month, indicates that we should use 
different parameters (monthly averages) to describe the con- 
ditions of heliosphere in the modulation process. In fact at 
100 AU, where particles are injected, the conditions of the 
solar activity are similar to those present at the Earth 14 
months before. Therefore, we can divide the heliosphere in 
14 regions as a function of the radius. For each region we 
evaluated fco, a and Vgw, in relation to the time spent by the 
solar wind to reach this region. We indicate the present treat- 
ment accounting for the time evolution of the solar parame- 
ters as a dynamic approach of the heliosphere. 



4 Results 

Results obtained with our propagation code are shown in 
Figs. |2]and[3] Simulated fluxes obtained using parameters 
dependent on the heliospheric region agree with measured 
data within the experimental error bars. This happens both 
in periods with A>0 (BESS-97), and in periods with A<0 
(PAMELA). This means that current treatment of the He- 
liosphere improves the understanding of the complex pro- 
cesses occurring inside the Solar Cavity. Our code can be 
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Fig. 2. Comparison of simulated - ratio at 1 AU and experimental 
data: BESS (1997). ^ 



10-^ 



,-e-io" 



10° 




#-f-*- 



LIS ratio 
Present Model 

PAMELA Adriani et. Al. 2010 



0.1 



1 10 
Energy (GeV) 



100 



^ http://omniweb.gsfc.nasa.gov/form/dxl . htrni] 
ihttp://www.sidc. oma.be/sunspot-data/, 



Fig. 3. Comparison of simulated - ratio at 1 AU and experimental 
data: PAMELA (2007-2008). ^ 

also used to predict CR fluxes for future measurements. The 
assumption is that diffusion coefficient, tilt angle and solar 
wind speed show a near-regular and almost periodic trend. 
The periodicity is two consecutive 1 1 -years solar cycles. We 
selected periods with a similar solar activity conditions and 
same solar field polarity of the time of interest; therefore ap- 
proximately 22 years before. We used the values measured 
in that periods as an estimation of the conditions of the helio- 
sphere. Simulations have been carried out in prevision of the 
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Fig. 4. Prediction of modulated ^ ratio at 1 AU for AMS-02. 



AMS-02 mission that will be installed on the ISS in 201 1 : we 
choose January 2012. For this period we show in Fig. |4]the 
predictions of GCR modulation for the p/p ratio. In order to 
reduce the uncertainty it is important to compare our model 
with the AMS-02 data because of the huge statistics and the 
long time covered. 



5 Conclusions 

We developed a 2D stochastic MC code for particles prop- 
agation across the heliosphere. We compared the ratios of 
p/p fluxes measured by BESS and PAMELA with those ob- 
tained from the present MC code. In the present calculations 
we used - for the parameters /cq, a and Vsw - values corre- 
sponding to the periods of data taking. This description of 
the heliosphere and the forward approach seem to properly 
account for the propa gation of GCR in th e solar cavity. Re- 
cent measurements l Adriani et al. (l2010l) 1 have pointed out 
the needs to reach a high level of accuracy in the modulation 
of the fluxes, in relation to the charge sign of the particles and 
the solar field polarity. This aspect will be even more crucial 
in the next generation of experiments like AMS-02. 
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